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Low and high intensity velocity selective coherent population trap- 
ping in a two-level system 
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Abstract. - An experimental investigation is made of sub-recoil cooling by velocity selective 
coherent population trapping in a two-level system in Sr. The experiment is carried out using the 
narrow linewidth intercombination line at 689 nm. Here, the ratio between the recoil shift and 
the linewidth is as high as 0.64. We show that, on top of a broader momentum profile, subrecoil 
features develop, whose amplitude is strongly dependent on the detuning from resonance. We 
attribute this structure to a velocity selective coherent population trapping mechanism. We also 
show that the population trapping phenomenon leads to complex momentum profiles in the case 
of highly saturated transitions, displaying a multitude of subrecoil features at integer multiples of 
the recoil momentum. 



Introduction. — The role of quantum interference in 
atomic physic has been demonstrated in several impres- 
sive experiments, e.g., in [1-3]. Already in the 1970'ies, 
it was shown that quantum interference can prevent ab- 
sorption in the presence of resonant light [4]. Indeed, if 
considering a A-shape configuration of internal states (two 
ground states and one excited state), there will exist a co- 
herent superposition of the two ground states, for which 
the two excitation amplitudes to the excited state interfere 
destructively. For counter propagating laser beams, such 
a dark state is velocity sensitive, and it can thus be used 
for cooling processes going below the single photon recoil 
energy. An experimental proof of such 'velocity selective 
coherent population trapping' (VSCPT) was first demon- 
strated in one dimension in 1988 [5]. The experimental 
signature consisted of two peaks in the momentum dis- 
tribution, centred at the momenta ±hk, and of sub-recoil 
widths. A few years later, VSCPT was also observed in 
2D and 3D [6,7]. Those experiments were done on the 
J g = 1 — * J c = 1 transition in metastable helium at 1.08 
/xm, using a near resonant laser with o~ + — o~ -polarisation 
configuration. Hence, two major requirements for VSCPT 
were fulfilled: firstly, the existence of a closed family of 
states with respect to the laser interaction; secondly the 
resulting dark state is also an eigenstate of the kinetic en- 
ergy Hamiltonian. Thus, the dark state gets an infinitely 
long lifetime leading to arbitrary narrow peaks in the mo- 



mentum distribution. In such conditions, very low tem- 
peratures have been obtained [8] . 

For transitions with higher angular momenta than in 
the above mentioned He*-experiments (i.e., with ground 
state angular momentum quantum numbers J g > 2), at 
least one dark state exists if (J = J g ) or (J c = J g — 1). 
However those dark states are not eigenstates of the ki- 
netic Hamiltonian. In reference [9], the authors suggest 
to add an extra off-resonant laser beam to induce light 
shifts in order to exactly compensate the kinetic energy 
mismatch. Without any compensation, the dark state ac- 
quires a finite lifetime due to motional coupling. As far 
as the lifetime of this state remains long with respect to 
other states, it can be favoured by the system [10]. Thus it 
may be possible to observe VSCPT with sub-recoil peaks 
at momenta zLMjhk with Mj > 1, where Mj is the mag- 
netic quantum number describing the projection of J g . So 
far no experimental data have been published using such 
a configuration. High momentum dark states has been re- 
ported for the case with a J s = 1 — > J c = 1 transition in 
metastable helium using a lin-angle-lin polarisation config- 
uration [11]. In that case, the dark states are characterised 
by two peaks at momenta ±Qhk, where Q is an integer. 
Thus, the high momentum states are still eigenstates of 
the kinetic Hamiltonian. 

From reference [10], we learn that observed momentum 
states are not necessarily totally dark states, and are thus 
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not eigenstates of the total atom-laser Hamiltonian. How- 
ever they should be the most protected states with respect 
to spontaneous emission processes. Following this idea 
and reference [11], the requirement of a closed family for 
VSCPT does not have to be strictly fulfilled. Indeed, con- 
sidering a two-level atom coupled with a nearly resonant 
laser, the coherent superposition of ground states with dif- 
ferent momenta; 

±(\ g; -hk)-\g ]+ hk)) (1) 

is not coupled to the excited state |e;0). However the 
three involved states, do not form a closed family because 
of the coupling to the states |e ; ±2hk). With a broad tran- 
sition, i.e., e = c<j r /r <C 1 (where T is the natural linewidth 
of the transition and uj r = hk 2 /2m is the recoil angular 
frequency), the coupling is not sensitive to a kinetic en- 
ergy mismatch between the |e ; 0) and |e ; ±2hk) states and 
VSCPT will not occur. If, however e w 1 (or even e >• 1), 
the unwanted transitions may be off-resonant, and the co- 
herent superposition in equation 1 acquire a long lifetime. 
This situation, leading to VSCPT, has been numerically 
studied in [12] and was observed on a helium beam, with 
e = 0.22 [13]. 

In this letter we report on observed VSCPT on the 
x So — > 3 Pi intercombination line in Sr at 689 nm. Here 
£ = 0.64 (with r = 2ir x 7.5 kHz). This more favorable 
value than the one used in [13] allows for a more pro- 
nounced observation of VSCPT in a two-level system at 
low intensity. The momentum distribution is still domi- 
nated by Doppler cooling effects, but an enhanced pop- 
ulation builds up at the momenta ±hk. The amplitudes 
of these peaks strongly depend on the atom-laser detun- 
ing. Moreover, for a higher saturation parameter, multiple 
peaks in the momentum distribution are observed. This 
structure is not linked to any closed family of states and, 
in contrast to the low saturation case, does not correspond 
to an eigenstate of the kinetic Hamiltonian. 

Experimental set-up. — The details of the experi- 
mental set-up, including two-stage cooling and trapping of 
Sr, can be found in [14,15]. Briefly, the strontium atoms 
are first accumulated and cooled in a ('blue') magneto- 
optical trap (MOT) on the (% -> x Pi) transition at 461 
nm, and are then transferred into another ('red') MOT, 
running on the ( 1 So — ► 3 Pi) intercombination line at 689 
nm, and with saturation intensity / S at = 3 ^W/cm 2 . The 
cold sample, with a 50 /im rms radius, contains about 10 6 
atoms at a temperature of 1 ^K. After the MOT phases, 
the atoms interact with an intensity balanced ID standing 
wave along a horizontal direction. A 0.1 mT vertical mag- 
netic bias field is applied in order to properly define the 
quantisation axis. This allows us to restrict the descrip- 
tion of the dynamics to two atomic levels. The polarisation 
of the counter propagating laser beams is adjusted to be 
parallel to the magnetic field axis. A typical duration of 
the ID optical molasses phase is between 0.5 ms and 2 



Fig. 1: Examples of time of flight images. The ID optical mo- 
lasses laser beams are along the horizontal axis in the figure, 
as well as in the experiment. From left to right, the angular 
detunings are S = — 3.1cj r , — 1.8w r , and — 0.4u; r , and the inten- 
sities are / = 57 aa t for all images, where 7 sat is the saturation 
intensity. 

ms. The interaction time is then always longer than the 
typical lifetime of the relevant internal states. Hence, the 
studies of the VSCPT cooling mechanism in this article 
are made in the steady state regime. However, this may 
not be the case for particular momentum distributions, 
where a steady state does not exist even for red frequency 
detuned laser [16]. 

The momentum distribution is extracted this distribu- 
tion using a time-of-flight (TOF) technique. The typical 
dark period of the TOF is 50 ms. Thereafter, a resonant 
probe at 461 nm is switched on for 40 /j,s, and the fluo- 
rescence signal is collected on an intensified CCD camera. 
The momentum resolution is limited by the finite size of 
the cloud after the ID optical molasses phase. If this phase 
is not too long, the typical rms radius of the cold cloud is 
still about 50 fim, which leads to a momentum resolution 
of hk/6. Images for background substraction are taken 
with the same procedure, but with no atoms in the blue 
and the red MOTs. 

Experimental results. — In figure 1, we show exam- 
ples of acquired TOF-images for three different detunings, 
but for otherwise identical laser parameters. The standing 
wave axis corresponds to the horizontal axis of the picture. 
By integrating along the vertical axis in fig. 1, we obtain 
more precise data for the momentum distributions. In fig- 
ure 2, we show such momentum distributions for a range 
of detunings. 

It is known that Doppler cooling on broad transitions 
(s <C 1) leads to Gaussian momentum distributions. In 
the present case, one can clearly see the more complex dis- 
tribution resulting from Doppler cooling on narrow tran- 
sitions. For instance we have observed that the mini- 
mum momentum dispersion is obtained at a detuning of 
about 5 = — 4w r , in contrast to the case of broad transi- 
tions, where the minimum dispersion occurs at <5 — — 0.5T, 
i.e., S ~ — 0.3w r (for e = 0.64) [16]. For a detuning of 
5 = — 0.5r, laser cooling on a narrow transition yields to 
two separated maxima, a 'double hump': the atoms are ex- 
pelled from the central region. This type of distribution is 
in qualitative agreement with the predicted non-stationary 
distribution reported in [16]. 

On the overall momentum distribution attributed to 
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Fig. 2: Momentum distributions for the intensity I — 5I sa t, 
and for angular detunings from +0.9 w r to — ll.lw r . These 
profiles are directly obtained from images as in fig. 1, by inte- 
grating the image density along the vertical axis, and plotting 
against the horizontal axis. The latter is converted to momen- 
tum units, taking into account the TOF expansion time. 
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Fig. 3: Momentum profiles from fig. 2, for the angular detun- 
ings (from top to bottom), — 4.4oj r , — 3.1lj t , and — 1.8a> r . The 
raw data is the red circles. Fits to four Gaussians (five in the 
lowermost case), are presented by a dashed black line. The 
blue triangles are the raw data, subtracted by the resulting 
fit function, with exception of the two Gaussians that appear 
close to ±p r . Thus, this is a good indication of the part of the 
atomic population that is in the semi-dark state. Vertical lines 
at ±p r are added, as a guide to the eye. 
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Fig. 4: Red squares: ratio of the atomic population in the 
VSCPT state, as a function of angular detuning. Blue line: 
ratio of the VSCPT feeding rate r Fccd and the VSCPT lifetime 

rVSCPT. 

Doppler cooling are superimposed two sub-recoil peaks at 
±Hk due to VSCPT [12]. The entire momentum distribu- 
tion is described as a number of Gaussians. Two of those, 
centered at p = ±p r , are used to account for VSCPT. In 
figure 3, we show such fits, corresponding to three of the 
traces in fig. 2. The figures also show the residuals from fits 
including only the broad Doppler features. The widths of 
the subrecoil VSCPT features are limited by the convolu- 
tion with the initial size of the cloud, and we can thus not 
extract relevant information about the VSCPT velocity 
distribution. However, by integrating the VSCPT-peaks, 
we get a measure of the fraction of the population that 
is in the long lived (VSCPT) state. In figure 4 we show 
the fraction of atoms in the VSCPT-state as a function of 
detuning. This quantity remains relatively small, within 
the few percent range, and peaks around — u> T below the 
atomic resonance. 

For higher laser intensity, the momentum distribution 
becomes more complex. An increasing number of subre- 
coil peaks, at momenta ±nhk (n > 1) are now observed, 
as shown in figure 5. As we will discuss in the next sec- 
tion, these subrecoil structures can also be attributed to 
a VSCPT mechanism. Those peaks are even less pro- 
nounced and broader than for the low intensity case and a 
quantitative description is not easy to achieve. It is how- 
ever important to note than the population fractions in 
the peaks are not necessarily the same or monotonously 
decreasing with momentum for a given experimental re- 
alisation. For example figure 6 clearly shows more pro- 
nounced peaks for the p = 0, ±2hk impulsion than for 
p = ±hk. We show in the following section that a long 
lived state appearing at high intensity can be attributed 
to the peaks at p = 0, ±2hk. 

Theoretical model and comparison with experi- 
ment. 

Effective Hamiltonian. A straightforward method to 
reveal the existence of long lived states is to diagonalize 
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Fig. 5: Time of flight images and momentum distributions for 
high intensity. From top to bottom, the intensities are / = 
30/ S at, 130 7 S at, and 350/ S at with 5 ~ — 5tv T . In the momentum 
distribution, vertical lines have been drawn at integer values of 

P/Pr- 



the effective Hamiltonian, which takes into account the 
finite lifetime of the excited state [12]. In the rotating- 
wave approximation, it takes the following form: 




Fig. 6: Time of flight image at / = 1307 sa t and S — —2uj r , and 
the corresponding momentum distribution. To improve the 
signal-to-noise ratio the profile is symmetrized with respect 
to the center of the distribution. The insert shows a resid- 
ual, where the slowly varying Doppler distribution has been 
removed. 
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Here il is the Rabi frequency of the atom laser coupling 
and we have ^ft- = -J— = sq , where I is the laser intensity 
and so the resonant saturation parameter. The coupling 
term in the Hamiltonian only connects ground and excited 
states that have momentum differences of hk. For this 
reason, one can reformulate the effective Hamiltonian in 
the following way: 



H c s — 
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(q + 2nhkf 
2m 
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|e ; q + nhk) (e ; q + nhk\ 
(q + 2nhkf 



2m 



\g;q + nhk) (g;q + nhk\ 



+ ~2 (\e;q + nhk) {g;q+ (n + l)hk\ 
+ |e ; q + nhk) (g;q+(n — l)hk\ ) ] , 



(3) 



with n being an integer and < q < hk. Each family of 
states is characterized by a q value and by an odd (resp. 
even) value of n for the ground state and an even (resp. 
odd) value of n for the excited state. 

Strictly speaking, each family contains an infinite 
number of members. However, one can remove high- 
momentum states, since the Doppler shift brings them far 



off-resonance. Hence to solve the eigenstate equation, we 
choose n max , a maximum value of n, such that the cou- 
pling between the |g(e) ; q ± (n max — l)hk) state and the 
|e(g) ; q ± n max hk) state is small and has any impact on 
the eigenstates of interest. In figs. 7a and 7b, we show 
the real and the imaginary parts of the eigenvalues, corre- 
sponding to some of the eigenstates of the effective Hamil- 
tonian (equation 3) for e = 0.64 with n max = 8 and q = 0, 
as a function of the Rabi frequency. 

Low intensity case. At low Rabi frequency, an ex- 
pected behaviour corresponding to VSCPT on a two level 
system is observed. The VSCPT state, namely the long 
lived state, has an eigenvalue which corresponds to the red 
line (full line in figure 7b). Indeed, its imaginary part re- 
mains small when SI increases, while its real part does not 
change significantly. At vanishing f2, we check that the 
eigenstate is, as expected, the one defined by expression 
(1). The smooth dependency of the VSCPT state lifetime 
is mainly due to the off-resonant coupling to the |e; 2±hk) 
states. Thus one has 
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T 2 + 4{5 - 3cj r ) 2 



(4) 



where r o g is the red dashed line in figure 7b. On the other 
hand, the green line (in bold in figure 7b) corresponds 
to the coupling state orthogonal to the VSCPT state for 
vanishing f2. At resonance, the lifetime Tq of this state is 
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Fig. 7: Real (a) and imaginary parts (solid lines) (b) of the 
eigenvalues of the effective Hamiltonian given by the expression 
3, with n max = 8 and q = 0, as a function of Q. Only the eigen- 
values corresponding to eigenstates for which the momentum 
\p\ < 2hk at SI = are shown. In (b) the dashed red and black 
lines correspond respectively to r o ff and T'^+uJt), whereas the 
dash-dotted line corresponds to Fc in units of F (see text). At 
17 ~ 4.5uj t and Q ps 12.5tJ r long lived states appear, pointed out 
by arrows and coloured circles. The histograms in (c) represent 
the populations of the relevant momentum states in the ground 
state (plain bars) and in the excited state (dashed bars). 



mainly due to the coupling to the |e; 0) state: 

m 2 



r — 

1 on — 



r 2 + 4(<5 + w r ) 2 



(5) 



Indeed r on , which corresponds to the dashed green line in 
figure 7b, coincides with the green line at small f2. 

One can estimate the population of the VSCPT state by 
a simple feed and loss mechanism. One then assumes that 
the population is proportional to the ratio Ff eec i/rvscPT, 
where Tf cc d ss (IVscpt + Tc)/2 is the feeding rate. Using 



the result given by the numerical simulation for evaluating 
rf ee d/rvscPT, we compare this ratio to the experimental 
data points in figure 4. The result is very satisfactory. 

High intensity case. If the two-level atom can be 
treated semi-classically, like for a broad transition for ex- 
ample, the two imaginary parts monotonously converge at 
high intensity to the same value, namely Im[£7] = — T/2. 
In the full quantum problem, similar behaviours are also 
observed. This is for example the case for the states 
\e;p = 0) and \g;p — 0) (black curves in figure 7b). In 
contrast, other states have unexpected non- monotonous 
behaviours. Some of them exhibit minima as a function 
of Q, sort as the one depicted by a green line at Q ~ 4.5a-> r 
and a blue one at Q w 12.5a-> r . Those states are of par- 
ticular interest as they are long lived ones. We will now 
demonstrate that those states are at the origin of the multi 
peak structures experimentally observed at high intensity. 

Figure 7c shows the normalized population of the mo- 
mentum states, for the two long lived state. If a two- 
peak structure at p = ±hk is a characteristic of VSCPT 
at low intensity, an increasing number of peaks is pre- 
dicted by our model and also observed in the experiment 
at higher intensity. For example the long lived state at 
fl « 4.5tc> r is mainly composed of three peaks, at p = 
and at p = ±2hk. At Q w 12.5w r a four-peak structure at 
p = ±hk and p = 3 ± Hk is expected. 

Lets focus on the simplest case, namely the three peak 
case at Q ss 4.5 w r , and derive an analytical expression for 
the state lifetime. In the view of the most abundant peak 
population, we restrict ourselves to n max = 2 and consider 
only the closed family {|g ; q) , |e ; q ± Hk) , |g, q ± 2hk}}. 
Under this condition, a straightforward diagonalisation of 
the effective Hamiltonian can be performed in the limit 
of high saturation. Here the kinetic terms and the radia- 
tive decay are removed and later been treated with the 
perturbation theory. One gets the following eigenvalues; 
A = 0, A = ±1/2, and A = ±V3/2 (in units of Hft), with 
the eigenstates 



|A = 0) 
|A = ±l/2) 
|A = ±V3/2) 



(\g;q-2hk)-\g;q) + 



1 

\g;q + 2Hk)) 
1 



(G) 



(-\g;q-2hk)T\e;q-hk)± 



± |e;<z + fifc) + |g;g + 2ftfc>) (7) 
' r \g;q~2hk) ±V3\e;q-hk) + 



2\g;q)±V3\e;q + hk) 
\g;q + 2Hk)) . 



(8) 



In this context the | A — 0) state has any excited state com- 
ponent with three peaks at p = and p = ± 2hk, thus very 
similar to the example shown in fig. 7c. The | A = 0) state 
gets a finite lifetime firstly due to the kinetic term which 
mix it to the |A = ±1/2) and |A = ±3/2) states and sec- 
ondly due to the off-resonant coupling to the |e ; q ± ihk) 
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states. Using perturbation theory, the imaginary part, of 
the |A = 0) state can be calculated. With q <C hk, one 
gets: 



(9) 



where 



Tkin = 

r o ff = 



(8m r ) 2 AT 

m 2 i 



i 



f-V 

\hk) 



3 2(5u> r ) 2 + fi 2 



1 + 



(2(5w r ) 2 + fl 2 ) 2 



\hk) 



(10) 



The general dependance of r^=o at q = 0, given by relation 
9, is fund to be in good agreement with the numerical sim- 
ulation presented in fig. 7. Moreover, the q dependence of 
Ta=o indicates that the long lived state is also velocity se- 
lective. This point is particularly important for a VSCPT 
cooling scheme. 

Momentum distributions corresponding to the |A = 0) 
state have also been observed in the experiment, as shown 
in fig. 6. The detailed VSCPT distribution is shown as a 
residual, where the slow varying envelop has been numeri- 
cally removed. We observe three peaks of the same height 
in good agreement with the prediction from fig. 7c. How- 
ever the measured intensity is I ~ 130/ sa t, i.e., i=a 12o; r . 
Even if the |A = 0) state is still a long-lived one at this 
value, the intensity is two times larger than the predicted 
optimum one (Q « 4.5w r ). This discrepancy may be due 
to an absence, in our simplified model, of a dynamical 
description of the population of the long-lived state. In- 
deed the model gives a prediction of the escape rate via 
the state lifetime but the feeding process is not described. 
The observation of those long lived states may occur at 
a larger or shifted range of intensity than the model pre- 
diction. Indeed only smooth changes, as function of the 
system parameter, on the multi peak structure have been 
observed so far. Moroever, since the image of fig. 6 does 
not reveal any peaks at ±3Hk, the ±hk structure can still 
be associated to the state (|g ; —fife) — |g ; +hk)) /V%, cor- 
responding to the low intensity case. 

Conclusions. — VSCPT-cooling on a two-level 
atomic system has been experimentally demonstrated. 
This is possible, since the atomic transition used, in Sr, 
has a narrow lincwidth, making the otherwise open fam- 
ily of momentum states, |g;— hk), |e;0), and \g;+hk), 
less open, due to kinetic detuning from other momentum 
states. Thus, a semi-dark state is formed. During the 
cooling, the long-lived state, of sub-recoil width, is fed 
with atoms, while there is simultaneously a constant loss. 
For a favourable ratio between the feeding rate and the 
loss rate, the momentum profile acquires narrow peaks, 
centered at ±hk, that lie on top of the Doppler cooled 
profile. At best, the steady-state population of the dark 
state reaches « 10%. 



At high saturation, complex momentum profiles arise 
with subrecoil peaks also resulting from a coherent pop- 
ulation trapping mechanism. With an analysis based on 
an effective Hamiltonian approach, we have identified the 
observed multipeak structures to long lived states. 
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